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Abstract 

We consider a system of reaction-diffusion equations with passive advection term and Lewis 
number Le not equal to one. Such systems are used to describe chemical reactions in a flow in 
a situation where temperature and material diffusivities are not equal, ft is expected that the 
fluid advection will distort the reaction front, increasing the area of reaction and thus speeding 
up the reaction process. While a variety of estimates on the influence of the flow on reaction 
are available for a single reaction-diffusion equation (corresponding to the case of Lewis number 
equal to one), the case of the system is largely open. We prove a general upper bound on the 
reaction rate in such systems in terms of the reaction rate for a single reaction-diffusion equation, 
showing that the long time average of reaction rate with Le ^ 1 does not exceed the Le = 1 case. 
Thus the upper estimates derived for Lc = f apply to the systems. Both front-like and compact 
initial data (hot blob) are considered. 

1 Introduction 

Systems of reaction-diffusion-advection equations describe numerous physical processes that occur 
when reactants not only diffuse but are also advected by a fluid or wind motion. The simplest 
situation is when the effect of temperature and concentration variations on the fluid motion may be 
neglected. Then evolution of temperature T{t,x) and concentration n(t,x) may be described by a 
system of two equations 

Tt+u-VT = KAT + ^g{T)n (1) 

rit + u ■ Vn = —An -g{T)n, x = (x, y), 

Le K 

where the velocity u is passive and presumed given. The nonlinearity g{T) is assumed to be of the 
KPP-type: 

g{0) = 0, g'{0) ^ 0, g{T) < g'{0)T for T > 0. 

The Lewis number Le, which is the ratio of temperature and material diffusivities may be different 
from one. The special case Le = 1 was studied in a great detail in the absence of advection beginning 
with the pioneering papers of Kolmogorov, Petrovski and Piskunov |15] and Fisher In this case 



for the initial data satisfying T+n = 1 the system (||) reduces to a single equation for the temperature 
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Tt + u-VT = KAT+'^fiT), (2) 



where f{T) = g(T){l — T). Recently the role of advection u in (^) has been a subject of active 
research (see [^] for a review, and |3|, @, |l^, 13, HI for very recent papers). In particular. 
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existence of traveling fronts in shear flows Q and pulsating traveling fronts in periodic flows was 
established |jl8|, ^. It was also shown that asymptotically front-like and compactly supported data 
propagate with the speed of the traveling fronts . The flow may have drastic effect on the front 
propagation, speeding up the reaction. The physical reason for this phenomenon is believed to be 
that the advection distorts the front, helping the hot material to warm up the cold one and increasing 
the area available for reaction. Various estimates for the speed of the traveling front in the presence 
of advection when Le = 1 were obtained in |^, ^ 14, 13 1, indicating the dependence of the 



burning enhancement on the intensity and geometry of the flow. 

An obvious restriction placed by the model (|2|) is the assumption Le = 1. There are some 
interesting situations where this is clearly not the case, such as a majority of reactions taking place 
in viscous liquids, or nuclear combustion in stars, where temperature diffusivity is much higher than 
material diffusivity and one can assume with a good degree of approximation that Le = oo. It is 
of interest, therefore, to compare the results already derived for (^) with what one can expect for 
the system case, where Le 7^ 1. The problem turns out to be much harder, in particular because 
of the lack of maximum principle for T. Surprisingly little is known about the system (|I|) when 
Le 7^ 1 even in the absence of advection. While it is known that traveling fronts exist, the set of 
allowed velocities may differ significantly from the single equation case, at least for some (non-KPP) 
reactions. A. Bonnet Q provided examples of reactions for which (P) has a unique traveling front, 
(H) has two disjoint intervals of possible traveling front velocities. Furthermore, it is not known that 
temperature T remains uniformly bounded in time. The best known estimate for ||T(t)||oo valid 
for u = grows like log log t for large t. The main purpose of this paper is to establish an upper 
bound on the bulk reaction rate for the full system (||). It shows that the rate of reaction for the 
system is bounded from above by the traveling front speed for the case Le = 1. Therefore Le 7^ 1 
does not provide a speed-up of reaction. This result is established both for front-like and compactly 
supported initial data. We notice that the papers ||7|, |l3| contain non-trivial upper bounds on the 
speed of front propagation in some cellular flows and in shear flows respectively. As a corollary of 
this paper, these bounds extend to the system case. 

For the front problem we consider the system (jl]) in a strip D = Mx x [0,H]y. The boundary 
conditions are periodic in y 

T{x,y + H)=T{x,y). (3) 
We assume that initially material on the left is burned, while on the right it is unburned: 

T{0,x,y) =To{x,y) ^ 1, n{0, x,y) = no{x,y) ^ as x ^ -00 (4) 

ro(x,y) 0, nQ{x,y) 1 as x — > +00. 

Initially temperature and concentration are equal to zero and one outside a finite interval: 

no(x, y) = 0, Tq{x, y) = 1 for x < —Lq and no(x, y) = 1, Tq{x, y) = for x > Lq (5) 

for some Lq > 0, and 

To{x, y) + no(x, y)<C, Tq > 0, 1 > no > (6) 
in [—Lq,Lq\. The flow u{x,y) G C^{1R'^) is periodic in x and y: 

u{x + L,y) = u{x,y), u{x,y + H) = u{x,y) (7) 

and incompressible: 

V • n = 0. 
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We assume in addition that it has mean zero: 

L H 



We define the reaction rate as 



dx I dyu{x, y) = 0. 

dxdy 



V{t)= / Tt{x,y)- 



H 



D 



and its time average as 



t 

{y)t = \j V{s)ds. (8) 



Long-time propagation speed is described by 

(y)oo =hmsup(y)t. 

t^oo 

In a way similar to |M] one may show that, because of the boundary conditions, we have 



V{t) = - I ntix,y)^ = ^ / g{T{t,x,y))n{t,x,y) 



dxdy 



D D 

It is known that if Le = 1 then there exist pulsating traveling waves of the form U{x — ct,x,y), 
periodic in the second two variables and monotonic in the first. Such solutions exist for c > c=k, and 
solutions with the initial data such as above propagate asymptotically with the minimal speed 
[10, 11]. Our main result for the front-like data is the following theorem. 

Theorem 1 Let c=k be the minimal speed of the pulsating traveling wave for Le = 1. Let T{t,x,y), 
n{t, X, y) he solution of ^) with periodic boundary conditions and front-like initial data, as in 
&' '^'^^^ arbitrary Lewis number Le. Then we have 

{y)oo < c*. 

Remark. In fact we prove that for any e > 0, there exists a finite constant such that for any t, 

{V)t<C,+£ + ^. 

Theorem |l| implies that difference in material and thermal diffusivities may not speed up the 
front propagation relative to the case Le = 1. 

Our second result applies to the situation when we have initially an isolated hot spot of material. 
Let To be compactly supported and uq satisfy @. We no longer impose periodic boundary conditions 
(^) but rather consider the problem in the whole space. Then the reaction rate is defined by 

Vit) = I T,{x,y)^. (9) 

Its time average is now scaled not as in @ but rather as 

t 
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and 

((y))oo = limsup((y))i, 

since the total area of an expanding blob of hot material grows as . We may also define a traveling 
front going in direction e solution of 

Ti + n • VT = kAT + —g(T){l - T) 

of the form U{e ■ x — ct, x) monotonic in the first variable and periodic in the second and such that 

lim C/(s,x) = l, lim U{s,x.) = 0. 

S—^ — O0 S— > + 00 

Such solutions exist for c > c*(0) with c*(0) being the minimal speed in direction 9. 

Theorem 2 Let T{t, x, y), n{t, x, y) he solution of ^ with the initial Tq that is compactly supported 
and uq that satisfies Then the long-time average reaction rate satisfies the upper bound 

{{V))oo<\j cl{9)d9. 



The first term in the above inequality is equal to the asymptotic bulk reaction rate of a solution of 
(|^) with To having compact support and Le = 1. Thus Theorem |^ is the analog of Theorem Q for 
such initial data. 

Finally, we note that our arguments also provide some information on the important case of 
ignition- type nonlinearity, where the function g{T) in (|l]) satisfies g(T) = for < T < Tq for some 
"ignition temperature" Tq < 1. Directly from the proofs, it is clear that in this case we get an upper 
bound on the reaction rate in terms of the minimal traveling wave velocity of the single equation 
d) with any KPP reaction f{T) satisfying f{T)>g{T){l-T). 

Theorem 3 Assume that the reaction function g(T) in is of ignition type. Let T{t,x,y), 
n{t,x,y) be solution of ^ with periodic boundary conditions and front-like initial data, as 
in ij^, with arbitrary Lewis number Le. Let fiT) be any reaction of KPP type (that is, 

positive on (0,1) and satisfying f'{0) > 0) such that f{T) > g{T){l — T). Let be the minimal 
speed of the traveling wave in with such f{T). Then we have 

{V)oo < C* 

Our methods also establish existence of a classical solution to (||) in a way similar to |^ extending 
the result of that paper to non-zero advection. 

We present the proof of Theorem || in the rest of the paper. The proof of Theorem |3| follows 
along the same argument. The main difficulty in the proof lies in the absence of known uniform L°° 
bounds on temperature T. If such bounds were available the proof would be greatly simplified. The 
proof of Theorem Q is very similar and hence is omitted. 

2 Outline of the proof of Theorem |I] 

The proof of Theorem |l| will proceed as follows. First we will show that the front cannot move to 
the right faster than c*: 
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Lemma 1 For any e > there exists a constant > and a constant > such that 

T{t,x,y) < c^e-^-(^-(^*+^)*). 

However, it is not yet known whether T(t, x, y) is uniformly bounded in time. Therefore we may not 
conclude from Lemma |l| that the conclusion of Theorem |l| holds since a lot of reaction may occur 
behind the front. To the best of our knowledge the strongest L°°-bound on T even for n = is 

r(t,x,y) <C(l + loglogt) 

obtained in It suffices for us to establish a weaker upper bound in the presence of non-zero 
advection. 

Lemma 2 For any e > there exists a constant > such that 

0<T{t,x,y) < Cec'K 

The proof of Lemma § essentially follows the ideas of Q with some modifications. Lemmas Q and § 
allow us to prove the analog of Lemma 1 for n(t, x, y): 

Lemma 3 For any e > there exists a function '^eH) ^ such that 

n{t, x,y)>l- ^s{x - (c* + e)t) 

and 



roc 

Jo 



Then Theorem |T] follows easily. 

Proof of Theorem [l] from Lemma ^. We first show that Lemma ^ implies Theorem ||. 
Observe that maximum principle implies that 



Then we have 



< n{t,x,y) < 1. 



1 f dxdy 
i^)t = 7 / X, y) - n{t, x, y)]- 



D 

(c,+e)t H oo H 

<- J dx J ^no{t,x,y) + - J dx J ^[1 - n{t,x,y)] 

-oo (c.+e)i 

c 

<c, + y (10) 
and Theorem |^ follows. □ 

3 Proof of Lemma |T| 

We now prove Lemma ^. Given any vector z G IR^ define the linear operator 

= k(V - z)2 - u • (V - z) + —g'(0) 

K 
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with periodic boundary conditions. The operator has a unique continuous positive eigenfunction 
(j){x,y) corresponding to a simple eigenvalue A(z). The positivity of the eigenfunction is standard, 
while setting i;^> = e"^, it is straightforward to verify that if u{x,y) is mean-zero and incompressible 

2 

then A(z) > -^g'{0). The first expressions for the propagation speed in terms of A(z) were given by 
Freidlin and Gartner |11, 10 1. The most convenient for our purpose representation for the minimal 
speed of propagation in direction e was given by Majda and Souganidis p^j: 



v{e) = inf 



A(ze) 



z>0 Z 



This expression also gives the asymptotic speed of propagation in direction e = (61,62) for Le = 1 
and general initial data satisfying (^), (|6|), (||). More precisely, we have for such data 



hm T(t, cet) 



t— >oo 



0. v{e) 

1, c < v{e) j ' 



ill) 



In particular the minimal speed of pulsating traveling front in direction ei = (1,0), that we denote 
by c^,, is given by 

A(zei) 



inf 



z>0 Z 

Given a number 2: > let ^{x, y; z) be the positive eigenfunction of Lz corresponding to A(2;ei): 



d_ 

dx 



+ 



9y2 



We define a comparison function 

cl,{t,x,y;z) = Cze-<''-'^^'^'/'^^{x,y;z) 
with the constant > so large that 

To{x,y) < (t>{0,x,y;z). 

We can find such C since ^{x,y; z) > ipQ > 0, and Tq vanishes for x > Lq. The function cj) satisfies 
a partial differential equation 

(f,t + u-V4> = KA(f> + —g'{0)(j) 

K 

Recall that n < 1, and g{T) < g'{0)T. Therefore T satisfies 

Tt + u-VT< kAT + —g'(0)T 

and so the maximum principle implies that for all t > we have 

T{t, X, y) < <t>{t, X, y; z) = ^e'^^^-^^^)*/^) vl;(;^, y- z). 

We choose then z > such that c*+e > A(z)/2; and obtain the conclusion of Lemma || since ^[x^y\z) 
is bounded from above. We remark that the proved bound remains true for any g(T) < MT with 
some M > not only KPP type. This observation will extend our arguments to give Theorem ^. □ 
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4 Proof of Lemma |2| 

Lemma ^ is proved using the technique of Q. We first prove a local L^-bound on T. 

Lemma 4 For any 7 > there exists a constant C > so that for any unit cube Q C we have 

TP{t,x,y)dxdy < Ce("^+^^')* 



with the constants a and (3 depending only on k and U = ||ii||oo- 

Proof of Lemma ^. Let 4>{t, x) be a smooth test function and F{T, n) be smooth in both variables. 
Then we have 



d_ 

dt 



(pit, x)F(r, n)(ix = / (ptFd^ 



(12) 



D 



D 



+ 



D 



FtkAT + ^FTg(T)n - Ftu ■ VT + ^FnAn - ^Fng(T)n - FnU ■ Vn 
K he K 



dx 



J[(t)t + u ■ V(p]Fd^ + ^ J {Ft - Fn)g{T)nd^ + ^ / ' 

D D D 



D D 

We rewrite the last expression above as follows: 



Ft AT + —An 
Le 



dx 



D 



Ft AT + —An 
Le 



(fx 



D 



FtVcp-VT +—V(P-Vn 
Le 



dx 



FttIVTP + FT„Vn • VT + — iVnP + — Vn • VT 

Le 



Le 



dx 



D 



and insert it into (|T^) to get 

j (pFdx = J [(l)t + u-V(j) + KA(f>] Fd:x: + k (^1 - J FnV(f) ■ Vndx 



d_ 

dt 



(13) 



D 



:^|Vn|2 + ("l + -L ) FT„Vn • VT + FttIVTI^ 



dx-^ f[Fn- FT]g{T)ndx. 

K J 
D 



We choose F in such a way that 
and 



Fn > 2Ft 



1 + ^) F^T<FnnFTT. 



(14) 
(15) 



Namely, we can take 
with A and e satisfying 



F = (^ + n + n'^)e^^ 
2e(A + 2)<l, ^>5(^1 + ^ 
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Then (^) and (15) hold and we obtain from (13): 

^ J (t)F<h^ < J + li • V0 + kAc/)] Fdyi + k{i- ^ J F„V(/) • Vndx 



(16) 



Le 



■|Vnp + FrT|Vr|^ 



d^-^ f Fng{T)nd^. 



D 



D 



We choose the function 6 of the form 



X 



(l+72|x-Xo|2)2' 



then 



|V(/>| < 27(/., \A(j)\ < 207^(/). 
Let us also denote U = ||n||oo, then we get from ([l^) with an appropriate constant C > 0: 



d 



dt 



— / (pFdx < C / C/7(/> + KY(t) Fdx + 2k 1 - — / Fnj(l)\Vn\dx - - 



Le 



K r .Fr, 



Le 



iVnPdx 



D 



D 



D 



D 



< 



(cU-i + ^7^) j ^Fdx + 2KLe -f^ J (py^dx. 



(17) 



(in the last step we replaced the quadratic expression involving |Vn| by its maximum). Observe that 

< 5e^^ < 2F 



Fl _(l + 2n)V^^^^^,^ 



2e=^ 



as long as ^ > 3. Then we have 



d_ 

dt 



J (f)Fdx <a J (t>Fdx 



(18) 



D 



D 



with 



Therefore we obtain 



with 



a = C 



C/7 + KJ 

j (t)Fdx < Ve."^ 

D 

(x)F(ro(x),no(x))dx< 



D 



Txe{A + 2) 
7^ 



Then we have for any unit cube Q and any positive integer k: 

cPFdx >Aj <pe^^dx >—j ^{x)T^dx > j T^dx. 

D D Q Q 

Here we have chosen xq in the definition of <j) to be the center of the cube Q. Therefore we obtain 

/ T''dx<A~^re''\l + j^)e^''kl 
Jq 



for positive integers k. Using interpolation we get for all p > 1 

Jq 

This finishes the proof of Lemma 

We prove now Lemma |2| using Lemma ^. Let G(t, x; z) be the Green's function of the advection- 
diffusion equation 

(pt + u ■ V0 = kA(/) 

posed in the whole space M^. Then we extend Tq periodically to the whole space from D and 
obtain 

t 

T(t,x) = J x; z)To(z)dz + J ds J G{t — s,x;z)g{T{s,z))n{s,z)dz 
IB? IB? 

t 

< I + J ds J G{t — s,x;z)g{T{s,z))n{s,z)dz. 



iR2 



The Green's function satisfies a uniform upper bound [|T^] 

G(t,x;z) < -^e-l"-"!'/^^'^*) =(5(t,x-z) 
with some positive constants C and v that depend on k and u{x,y). Then we get 



(19) 



r(t,x) < 1 + 5(0) y y G(t- s,x-z)T(s,z)n(s,z)(iz 

iR2 



1 + 



gg^(O) 



j ds j e-l^~^l'/(^^(*-^»r(s,z)n(s,z)(iz. 



JR2 



We split the last integral into sum of integrals over unit cubes Qj with the cube Qq centered at the 
point x: 



j e-l'^-^l'/(Mt-^))r(s,z)n(s,z)dz = ^ j e 
IB? ^ Qi 



|x-z|2/(4Kt-s)) 



T(s, z)n(s, z)dz 



-{dist{x,Q,))V(8i'{t-s)) /■ |x-z|2/(8;^{t-^)) 



^T(s, z)n(s, z)(iz 



< 



in 


the integral above: 














< 


J g-q|x-z|2/(8;.(t-^))^2 




y TP{s,z)nP{s,z)dz 




- 7R2 




Qj 



l/p 



87rz/(t - s) 



1/-? 



TP{s,z)dz 



i/p 



< C{j,i^,q)Leit- sy/'^e 



l/qja'y+l3-y^)t/p 



Therefore we have 

G(t — s, X — z)T{s, z)n{s, z)dz < 



C 



Note that if z € Qj with j 7^ then 



and therefore 



{A'Kv{t - S)) 



[t - s) V'?e(a7W)*/p ^ g-(dist(x,Q,)V(8j.(t-s)) 



|x — z| < 2dist(x, Qj) 



G{t - s, X - z)T{s, z)n(s, z)dz < C{t - sY'i'^e 



IB? 



1 + j e-l"-^l'/(^6.(t-s))^^ 

iR2 



< C{t - s)-i/Pe("^+'^^')*/P [1 + C{t - s)] < C7e("^+^^')*/^' [{t - sY^'p + {t - sf/'' 
Finally integrating over s £ [0, t] we obtain 



r(t,x) < 1 + Ce^^^+'^^'^^/P 



^l/g , ^1+1/, 



However, 7 > is arbitrary and thus Lemma |2| follows. □ 

5 Proof of Lemma 3. 

Let us define W = 1 — C . It satisfies a differential equation 

2 

Wt + u- VW = ^AW + —g(T)n. 

Le K 

Lemmas || and |2| imply that given any e > we may find Cs, Xe > such that 

T{t, X, y) < Cee'\ T{t, x, y) < C,e-^^("-^^*) 
with = c=K + e. Let us define 

2 

Rit,x) = ^min(C7,e"*,C,e-^^(^-"=*)) . 

Then we have by the maximum principle 

W{t,x,y) < ^{t,x,y) 
with the function $ satisfying the initial value problem 

$f + li • V^> = -^A^ + R(t, x) 
Le 

$(0, X, y) = To{x, y) < H{-x + Lq). 
Here H{x) is the Heaviside function: 

0,x < 



H{x) 



l,x > 



10 



and Lq is as in (^). In the sequel we will assume without loss of generality that Lq = 0. We let 
r(t,x;x') be the fundamental solution of 



Le 

and obtain (x' = {x',y')) 



$(t,x) = j r(t,x;x')F(-2;')(ix + j ds j r(t - s, x; x')i?(s, a;') = + ^2 



iR2 iR2 



We will bound and $2 separately. We have a bound for r(t,x;x') similar to (19): 

r(t,x;x') < — e-l''-'''l'/(^''*) (20) 

Kt 

with the constant C depending on k, u and Le. Then we obtain 

—x/\/ Cnt 

for 2; > 0. Observe that given any C > we may choose Aq > and B > so that for all x > and 
all t > we have 

Indeed we need to find Aq such that 

n2 ~1 



min < Xqx + XoCet )■ > —00. 

x>0,t>0 I CKt I 



However, we have for x > 0, t > 



X 



2 



2x^yXoCe 



Xqx + XQCet > ""^ '^""^ - Aox > 



for Aq sufficiently small. Therefore ( |2l| ) holds for such Aq and thus we have 

^i(i,x) < Ce-^«(^-^^*) for x > 0. (22) 

Our next goal is to obtain such bound for $2. We use the inequality (20) to get 

t ^ 

^2(t,x) <C fds f _4L=e-("-"')'/(^'^(*-^))i?(s,x'). (23) 
J J \ n\t — s) 
M 



Recall that R{t, x) is defined by 

R{t,x) 

Here X{t) is defined by 



Cee^\ x<X{t) 

(J^^-\4x-Cet) ^ X > X{t) 



et = -Xe{X{t)-Cet), X{t)=[c,-j-]t. (24) 
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We split the integral in (|23|) accordingly: 

t X{s) 

^/i^{t - s) 



Mt,^)<C',fds I , -{x-x'Y/{CK{t-s))^es (25) 

J J \l K,\t — S) 

t oo 



-oo 
X{s) 

The term $21 is bounded as follows: 



t {X{s)-x)/^CK(t-s) 

$2i(t,x) <Ce j ds j dx'e-^"e" 



-00 

2 , , „,2 



Recall that for a < we have /"^^ e ^ dx < e " . Therefore we have for x > Cet: 

t t 
$2i(i,x) <C, [ dse"-(^W-^)'/(^'^(*-*)) <C, [ dse"-(^-'==*)'/(^'^(*~^)). 



We use now the inequality (^T]) to get 

t 

$2i(t,x) <CeJ dse"-^o(^-'=^'-'=='(*-')) < Ce-^o^^-^^^+'/^o)*). (26) 


It remains to bound <I>22- The function <I>22 solves the initial value problem 



dt dx"^ 

«'22(0,X) =0 



+ C2e~^^^^~'^'^H{x - X{t)) 



with H{x) being the Heaviside function. However, we have using (24) 

e-^e{x-c,t)jj^^ _ < e''H{-x + Qt) + e-^^(^-"^*)/7(x - c,t), 

and therefore 

^•22 < $23 + $24. 

The function $23 solves the initial value problem 

$23(0, x) =0. 

It is bounded exactly in the same way as $21 but with X{t) replaced by x{t) = Cet. This gives an 
upper bound 

$23(t,x) < C7^e-^°(^-(^-+^/^o)*). (27) 
The function $24 solves the initial value problem 



dt dx"^ 
«>24(0,a;) = 



12 



It may be bounded from above by 

^2dt,x) < ^{x-Cst). 

Here is a positive solution of 

A general solution of the above ODE is given by 



Ck 



1 _ e-c.€/(C«) 



for ^ < 



and for ^ > 

Ck 

^(0 = ^0 + *i 1-e 



Let us require that 



Then the function ^ takes the form 



+ 



Co 



+ 



Co 



XeCK 



1-e 



Cn , C2 
^0 + ^1 + ^ = 0. 



for ^ < 0, and 



^(0 



Co 



XJcr - CnXr, 



-Cei/iCK) 



+ 



Co 



Xe{Ce - CkXs) 



for ,^ > (we can always assume that the denominator is nonzero by shifting a little bit). Therefore 
^' > as long as ^0 is non-negative, in particular we may choose ^'o = 0. Then we have 



00 

^24{t, X) < ^{X - Cet), J ^{Od^ < 



00. 



(28) 



This, together with (22), ( pq ) and ( p7| ) finishes the proof of Lemma ^.n 
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